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16, Abstract 

The amplification or reduction of unsteady velocity perturba- 
tions under the influence of strong flow acceleration or de.- 
celeration is studied. Supersonic flows with large velocity 
and pressure gradients are examined and the conditions in 
which the velocity fluctuations depend on the action of the 
average gradients of pressure and velocity, rather than turbu' 
lence, are described. Experimental results obtained on a 
cylinder positioned paraj.lel to the flow are analyzed statis- 
tically and interpreted as a "return to laminar" process. It 
is shown that this "return to laminar" implies negative value; 
in the turbulence production terms for kinetic energy. A 
simple geometrical representation of the Reynolds stress pro- 
duction is given. It is also shown that in the experimental 
case the influence of viscous energy dissipation cannot be 
neglected. 
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THE INFLUENCE OP A HIGH PRESSURE GRADIENT ON 
UNSTEADY VELOCITY PERTURBATIONS IN THE CASE OF 
A TURBULENT SUPERSONIC FLOW 

J. P. Dussauge and J. F. Debieve 


I. Introduction /I* 

The amplification or reduction of stationary velocity 
perturbations under the effect of strong accelerations or decele- 
rations of a permanent average flow is studied. The cases examined 
relate to supersonic flows in which there are high velocity and 
pressure gradients. In these conditions, the evolution of the 
turbulent velocity fluctuations processed by the statistical 
method, presents special aspects, which we attempt to analyze 
here. First of all we will define the conditions in which the 
evolutions considered are not entirely specific of the turbulence, 
then the experimental results obtained in a flow undergoing 
expansion are described. These results, which have already been 
published are analyzed under a new aspect: they are compared with 

a relaminarization criterion; the latter is interpreted in terms 
of "production of turbulence", of which a few analytical elements 
are given. Finally we make an estimate showing that in the 
studied expansion, the effect of the non linear and viscous terms 
involved in the kinetic energy equation of turbulence does not 
appear predominant. 


II. Analysis of the Problem 

We consider a permanent average flow, on which unstationary 
fluctuations of temperature, pressure, volumic mass and speed are 
superimposed. The le’''el of these fluctuations is low as compared 
with that of the average quantities; furthermore it may vary under 
the effect of gradients of average velocity and pressure. In this 

*Numbers in the margin indicate pagination in the foreign text. 
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study, we are particularly interested in the velocity fluctuations, 
which will be processed statistically, and in their evolution as 
a function of the variations of the average field. In certain 
conditions, this evolution may show aspects which are not totally 
specific of turbulence. The object of this paragraph is to 
specify these conditions. 


Let us consider the equation of kinetic energy of turbulence, /2 
established with the use of average mass weighted quantities. 
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represents the space variables, the velocity components, f 
the volumic mass, /Jv the pressure. The symbol represents an 
average, a mass weighted average and ' a fluctuation. The mass 
weighted average of a quantity is defined by 

p OO- 'UJ' which implies z O • 

The different terms of equation (1) may be interpreted as 
follows (J. Gaviglio, J. P. Dussauge, J. F. Debieve, A. Favre, 

1977) : the terra I represents the variation of along an 

average current line; II is a terra of "kinetic" production; III 
represents an "enthalpic" production; the terra IV, for a solenoidal 
velocity fluctuation field ( 9 - O) is usually grouped with 

the diffusion terra V; the terra VI designates energy dissipation 
under the effect of viscosity. 

Seeing that in the accelerated or decelerated areas in 
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supersonic flow, p may vary greatly, we will consider the evolution 
of pu.',a'. j\r , which is equal to w/.; if the term ^Ui;, uu'7 j p* of /3 

the third order with respect to the fluctuations, is negligible. 

In the turbulent flows which undergo little change, diffusion 
and dissipation play an important role. For example, in a boundary 
layer in equilibrium, the terms of production, diffusion and 
dissipation are of the same order of magnitude. The latter is equal 
to , if it is assumed that the dissipation has the same 

form in compressible and incompressible flows, q* ~ and 

-A is a space scale characteristic of energy bearing turbulent 
structures, for example a space scale derived from correlation 
(Tennekes and J. L. Lumley, 1972). 

In a distortion region, that is, wherever the average flow 
varies to a considerable extent, the order of magnitude of the 
production terms may be changed. It happens that for moderate 
Mach numbers ( M ^*^), the terms of kinetic and enthalpic production 
have comparable values (J. P. Dussauge, J. Gaviglio, A. Favre, 

1978) . Their order of magnitude may be estimated by p /Lei , 

in which and Lj designate a characterxstic scale of the evolution 
of the average velocity and a space scale relating to distortion 
respectively. 


If the terms of diffusion and dissipation do not vary much 
in the distortion, they can be disregarded, on the condition that: 
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which is equivalent to writing 
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We may estimate the condition for which the dissipation rate will 
not vary much: it is sufficient that the time during which a 


fluid particle undergoes distortion should be less than its 
"characteristitii time scale” (or "lifetime”) o^f-^ • Therefore 
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(see for example, Batr-aelor, 1967) 


U is an estimate of the local average velocity within the 
distortion. As regards the diffusion term, it must be considered 
separately in each special case (for example, see Paragraph V. 1.). 

If the inequalities (2) and (3) apply and if the diffusion 
rate varies little, this would mean that the flow is such that 
the influence of the effects of the diffusion and dissipation terms, 
specific of turbulence, is negligible. In this case, we find a 
type of flow in which the evolution of the level of unstationary 
fluctuations is relatively independent of the classical properties 
of turbulence. 


III. Studied Flow 

III. 1. Description 

The flow in which the variations of the statistical level of 
the fluctuations were observed is the close wake developing 
downstream of a truncated cylinder placed parallel to the flow 
(Fig. 1) . This configuration was already described in previous 
publications (J. Gaviglio, J. P. Dussauge, J. F. Debieve, A. Favre, 

1977, J. P. Dussauge, J. Gaviglio, A. Favre, 1978). Its main 
characteristics are recalled below. 

The pressure generating the potential flow is 0.375 atmosphere. /5 
The Mach number is 2.3 in the external flow, upstream of the base 
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of the cylinder. A fully developed turbulent boundary layer 
develops orv the cylinder. The boundary layer, near the end of 
the cylinder, is subject to expansion. There is detachment, 
with the constitution of recirculation area surrounded by a 
mixing zone. Recompression takes place during which there is 
confluence of the mixing area. The wake is formed downstream. 

The region which we will be discussing chiefly is the 
boundary layer undergoing expansion. 

III. 2. Measurements Carried Out 

III. 2. 1. The results described are of two types: 
measurements of average quantities by pressure probes and measure- 
ments of standard deviations and correlation coefficients relating 
to velocity and temperature fluctuations. In the following, the 
equations are written in a system of coordinates bound with the 
average current lines. The direction s is taken along their 
tangent, y along the normal, </' along the binormal. The only 
component of the avc.rage velocity is tv and the velocity fluctuations 
in the three previous directions are noted u' , 0 

designates the temperature; .* the distance in a given point and 
the axis of the cylinder of diameter cl ; ^ refers to the average 
current line of zero rate of flow. 

III. 2. 2. Initial Boundary Layer 

The characteristics of the boundary layer upstream of 
the expansion are as follows: its thickness S(o) is 4 mm; the 

Reynolds number referred to the thickness of the quantity of 
motion is 10^. The friction coefficient, derived from the average 
measurements is ~ 2,9. io“^. % is the friction 

on the wall; the index tsa* relates to the external flow. 

Standard deviation measurements of velocity and /6 

temperature fluctuations were carried out with a hot wire anemometer, 
of the "constant intensity" type. Because of difficult measurement 
conditions, the measurements were corrected by an adapted method 
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(J. GavigXio, J. P. Dussauge, 1977) and compared with the known 
results relating to boundary layers in equilibrium (Fig. 2) . 

The measurements carried out do indeed correspond to conditions 
of boundary layer in equilibrium. 

III. 2. 3. Expansion Area 


In the expansion area, measurements of pressure and, 
average velocity were carried out, making it possible to estimate 
the gradients of these quantities. An example is given by 
Figure 3, which presents profiles of average velocity. In the 
expansion, the results of the measurements were compared with 
the calculated profile assuming that along an average current line, 
the evolution is isentropic and the friction forces are small as 
compared with the pressure forces. In these conditions the first 
equation of the quantity of motion is written, in a system of 
coordinates bound with the average current lines: 


pa 




It is apparent that there is satisfactory consistency 
between measurement and calculation. The shape of the profiles 
obtained shows the effect of the expansion on the average velocity 
gradients: in particular there is creation of a positive gradient 

d’^l'do and considerable reduction of 

As regards the pressure gradients, we find a negative 
longitudinal gradient Qpybs and a positive transversal 

gradient. 

Measurements of standard deviations of velocity 
fluctuations were carried out in expansion (Fig. 4) . The results 
show that in this area there is a considerable drop in the level 
of the longitudinal velocity fluctuation. Analytical elements 
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which make it possible to relate this v^volution with that of the 
average gradients are proposed in paragraphs IV, 1. and V. 


III. 2. 4. Results Common to Any Plow 
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The measurements of temperature fluctuation and of 
the correlation coefficient between the longitudinal component of 
the velocity and temperature show that in any flow, including 
accelerations and decelerations; 

1. the formula derived from "Strong Reynolds Analogy" 
applies (Pig. 5) : 






M is the Mach number and V* the ratio of the specific heats of 
the gas . 

2. The correlation coefficient between velocity and temperature 
is practically constant and its value is close to -0,8 (Pig. 6) 
which is close to the values measured in heated subsonic flow. 

These results will be used in paragraph IV. 1. 2. 


IV. Interpretation of the Results of Measurements 

* 

IV. 1. Reninders 


IV. 1. 1. Criterion of Narasimha and Visvanath 

Since the studies of Sternberg (1954) , we know that 
boundary layers exposed to considerable accelerations, especially 
in centered expansions, in supersonic flow, have a tendency to 
"relaminarization. " This relaminarization is expressed in a 
decrease of the intensities of turbulence, and a variation of the 
value of the parameters like the parietal thermal factor. Narasimha 
and Visvanath (1975) reviewed the experiments in which relaminari- 
zation could occur, in view of establishing simple criteria making 
it possible to foresee the cases of "relaminarization." 
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These authors have reasoned as follows: they assume that in the 

expansion, the friction forces are small as compared with the 
pressure forces. Purthernrare, they take i^'to account the fact £8 

that in the expansions observed, the length over which the latter 
are spread are of the order of the thickness ^ of the boundary 
layer. The longitudinal gradient may therefore be estimated as 
-A'p /£ , in which ip the absolute value of the pressure 

deviation on either side of the expansion. If the pressure forces 
are great enough, relaminarization takes place. The ratio of these 
two forces may be estimated as follows: 

-to/s 

'to is the wall friction of the initial layer. 

Reviewing the measurements conducted by different 
authors, Narasimha and Visvanath find that for , 

there is return to the laminar state, and that for £jv/T;o<50 • 
this return is not mentioned. 


IV. 1. 2. Production of Turbulence 


The decrease of velocity fluctuations in the expansion 
of the near wake flow has been attributed in previous publications 
(J. Gaviglio, J. P. Dussauge, J. P. Debieve, A. Favre, 1977; J. P. 
Dussauge, J. Gaviglio, A. Favre, 1978) to negative values of terms 
of "production of turbulence", present in the equations of kinetic 
energy of the turbulence and in the equations of the Reynolds 
tensions. Let us recall briefly the main results. 


The terms of production are given as: 



J 
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* 

P will be designated as the simple term of "production'*, 
to differentiate it from kinetic and enthalpic productions. 

P its estimated with the following hypotheses; 

a. The field of velocity fluctuations is 
solenoidal /dzi^O) . 

The temperature fluctuations are practically 

(Laufer 1969) 

The values of the correlation coefficients 
between the standard deviations of the different velocity components, 
between the fluctuations of velocity and temperature are not very 
different from what they are in the boundary layer in equilibrium. 

d. The values of the ratios between the 
standard deviations of the fluctuations of the different velocity 
components are close to those found in the equilibrium boundary 
layers. 

e. The standard deviations of the longitudinal 
fluctuations of velocity and the temperature fluctuations satisfy 
the relation given by the "Strong Reynolds Analogy" mentioned in 
paragraph III. 2. 

Furthermore the terms of production of turbulence 
relating to each of the variances of velocity fluctuation 
*uy<-, , have been given in a coordinate system bound with 

average current lines. Taking into consideration the average 
geometry of the flow and the direction of mass transfers by 
turbulence, it may thus be shown that the absolute value of the 
production of is negligible as compared with those of 

, and ; furthermore, the production of is always 

negative. Therefore in this case if the production of 
is negative, the same will hold good for the production of 

• (Dussauge, Gaviglio, Favre, 1978.) 

The production of turbulence relating to 

is written ; 


b. 

isobar ic; , 
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The only positive term is " y^, But Figure 2 shows that in 

the expansion, ’S'^lon decreases considerably. With the hypotheses 
put forward previously, it is found that the terms diZ /d^ and 

are of opposite signs and of the same order of magnitude /lO 
in absolute value. The negative term Pu? is then highest in 

absolute value, it gives than its sign to ^ . Furthermore we 
found consistency in the order of magnitude between the estimate 
of V. and the measured rate of decrease of *.1'* . Thus the 

evolution observed is explained to a considerable extent by the 
negative values of the term of production of turbulence. 


IV. 2. Interpretation of the Criterion of Relaminarization 
of Narasimha and Visvanath as a Function of the 
Production of Turbulence 


The guiding principle of Narasimha and Visvanath in the 
establishment of their criterion is that the pressure forces are 
large as compared with friction forces caused by turbulence, in 
an expansion. This is an argument presented in terms of quantity 
of motion, which provides little information on the evolution of 
the kinetic energy of agitation. It is shown here that the 
criterion of Narasimha and Visvanath may be associated with 
negative values of the production of > therefore production of 
turbulence. Let us consider again the expression of relating 

to . 

P 

p- 

* ” F’ “Dg dT~ 


If the friction forces are neligible as compared with the 
pressure forces, we have (Paragraph III. 2. 2.); 
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Using the relations: 


£ = ^ C H*- 

Q .'t Vi 


and 
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The second term is negative. is also negative if: 
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Because of the geometry of the expansion area, we have 

, in which is the difference of pressure on either 

side of the expansion. 


Likewise 
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It was noted that is probably an upper bound of _ ^uT^P’ , and 

upper bound of dit/0^ (Fig. 3), outside the wall 

law: 

^ 

u.'*' 



We find: 


For , r = -O.D and \|u’^ - d.'IO*^' 

V4r *• 


£± > CO 
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If on the other hand, we calculate , in the case of 

our experiment, we find the value 46; we revealed here a negative 
production of turbulence. The condition appears /12 

therefore to be superfluous. This is no doubt due to the 
considerable excess taken when it is assumed that is of 

the order ^/S . Actually Fig. 3 shows that in the accelerated 
region, we would be having rather O,^ Uv /j . Furthermore 

the value 70 given by Narasimha and Visvanath is really of the 
same order of magnitude as the two above- indicated estimates, 
which correspond to negative values of the production term. 


IV. 3. Remarks on the Term of Kinetic Production 


In the preceding, we discussed mainly the production of the 
kinetic energy of turbulence. But the analysis of the behavior of 
this term led to considering the production of the normal Reynolds 
tension pu'* . A few characteristics are given here of the effect 
of the kinetic production of Reynolds tensions. It is found that 
it is difficult to evaluate the latter. But it is possible to 
obtain a simple geometrical representation of it taking into 
account the geometry of the average flow. 


To do this, let us- consider the quantity Q, r 
representative of the agitation in the direction of a vector .■ 
The evolution with the average movement according to the law 
(Debieve, 1978) is imposed on the vector tJ ; 


^ f.'; = 






with 


- i- 
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In the case of a pure deformation, this corresponds to the 
ordinary concept of the transport of a vector by the average 
movement: would be carried away by the fluid with contraction 

or dilation according to the deformation of the average raoveroent. 
The evolution of Q is given by; 
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It was shown (J. F. Debieve, 1978) that we can consider only the 
part of the tensions whose evolution is caused by production, 
specifically: 









u. 


= - 7^, f .1. u. 
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We obtain therefore: 



is a constant 


fu'n /f corresponds to velocity fluctuations in the direction . 

Thus we obtain a decrease of the level of the velocity 
fluctuations in the directions which dilate and conversely. As 
an example, let us examine the case of a single longitudinal 
gradient of stationary velocity (expansion or compression) . In 
the transversal directions, the vector remains constant, the 
transversal velocity fluctuations are unchanged. On the other 
hand, the longitudinal fluctuations will increase or decrease in 
the ratio We will therefore observe a decrease of 



in an expansion and an increase in compression. 


These results reveal simply how the different Reynolds tensions 
may vary one with respect to the other under the effect of spatial 
variations of the average velocity field. 


V. Discussion of the Validity of the Approximations Taken 
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V. 1. In paragraph II we showed ihe conditions in which it 
was possible to disregard the specif » tfects of turbulence. Two 
conditions were given: 

^ y. 

AU A.. 




H «4 

U A-y 


The first inequality provides the assurance that we may 
disregard the non linear terms and those expressing the viscous 
effects, with the reservation that they retain the same order of 
magnitude in distortion, and upstream of the latter. If the 
second inequality applies it may be estimated that the dissipation 
rate will not be changed in the distorted area. 

B f> q ' fc tTT 

The order of magnitude of the diffusion term I ‘ 
may be estimated: 


0 

dt-i 



» 





If the expansion considered were powerful enough to suppress all 
fluctuations, we would have: 

_ 

which probably furnishes an order of magnitude of the term 
considered. 
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it may /IS 


On the other hand, as regards the term -..j 
be considered that it represents a diffusion of the kinetic 
energy of the turbulence by the velocity fluctuation v'. 

This diffusion takes place over a distance comparable to an 
average dimension of the turbulent perturbations considered, for 
instance, the integral scale derived from spatial correlations or 
autocorrelation . We have therefore 






\j* 


A* 




In the studied case we find; Lj, tT . if, moreover, 
varies little in the expansion, we have Lj • The 

two diffusion terms are therefore really of the same order of 
magnitude as upstream of the expansion. 


V. 2. The application of the previously described inequalities 
implies that one should know a characteristic scale of the energy 
carrying turbulent perturbations _A.^ . 

There are many measurements of these scales in subsonic flow, 
but little information is to be found in the literature on these 
scales in the supersonic boundary layer. Demetriades measured 
them in the case of a wake, with Mach number 3, and in a cooled 
boundary layer, Mach number 9. Preliminary measurements were 
therefore carried out at the l.M.S.T. (Institute of Statistical 
Mechanics of Turbulence) in the case of supersonic boundary layers. 


The determination is carried out on the basis of the measurement 
with hot wire anemometer of the autocorrelation coefficient • 

V7e find by integration an integral time scale 

_A. . :: 

t. / 

J\.^ corresponds to the time of passage of a perturbation of 
magnitude -A., 4.. If it is assumed that these perturbations undergo 
convection through the average flow with the local average value 
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vC we find 


A^= K 


In the near wake flow, autocorrelation measureinenti! of the /16 
signal given by the hot wire anemometer were carried out for an 
overheating coefficient 0.2. In this case, the anemometer is 
sensitive both to speed and temperature. The integral scales 
found are of the order 0.1cf(J. Gaviglio, J. P. Dussauge, J. F. 

Debieve, A. Pavre, 1977) . 


More detailed measurements have been conducted in a boundary 
layer on plane plate for Mach number 1.8, by applying to the 
correlation measurements the technique of the "fluctuations 
diagram", which permits the separation of che respective effects 
of velocity and temperature. The preliminary results indicate 
that the integral velocity scale measured should be between 0.2 S 
and 0.4 (if . We recall that in subsonic boundary layer in equilibrium -A. 
is 0.45 (S' . With the results of the first measurements, we find 

that the parameter q* Ud is between 0.1 and 0.3. It appears 

U Am, 


therefore that the dissipation rate hardly changes in the expansion 
area. If it is considered that the scale is of the order of 
the increase of the longitudinal velocity component in expansion 
~~ "yUK order of 1. Actually if we carry out a more 

detailed estimate of the term of productj.on in expan sio n at a 
distance (r,. r,Vf ~ we find that “ ■** f*—/ 

is of the order /CCo) • The drssxpation term assumed to 

be equal to the production terra in the initial layer is of the 
order 0.1 . 


It seems therefore that in the expansion, the level reduction 
due to dissipation is not negligible. But it appears that the 
reduction observed in the levels does seem to have been caused 
mainly by the production of turbulence. 


Another case may be examined; that of the compression undergone 
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by the mixture layer near the reassembly (Fig. 1). This compression 
is spread over a fairly long distance (2 to 3 times the thickness 
of the turbulent layer) and we find in it considerable velocity 
fluctuations whose intensity is of the order of 15%. In this /17 

case the time of passage to compression is of the same order of 
magnitude as the characteristic time of large turbulent structures. 
Therefore the diffusion and dissipation effects must be taken into 
consideration to describe even qualitatively the evolution of the 
different turbulent flows. 


VI. Conclusions 

On the basis of the considerations on the characteristics 
of the turbulence of clipped flows in the presence of high gradients 
of average velocity and pressure, it was possible to reveal 
conditions in which the evolution of the velocity fluctuations 
does not depend totally on the specific properties of turbulence, 
such as diffusion and dissipation, but especially on the effect of 
the gradients of the Average quantities. In an adapted experimental 
case, the evolution of the statistical level of the velocity 
fluctuations is interpreted as a function of the gradients of the 
different average quantities. It was possible to relate this 
analysis to a relaminarization coefficient. It was shown in 
particular that a "relaminarization" is probably accompanied by 
negative values of the term of turbulence production of kinetic 
energy; a simple geometric representation of the effect of the 
production of Reynolds tensions in the average velocity gradients 
is given. Finally it was shown that in the case of the studied 
experiment, although the effect of the production of turbulence 
is predominant, the effect of energy dissipation by viscosity should 
not be disregarded. 

Institute of Statistical Mechanics of Turbulence 
12, avenue Gdndral Leclerc 
13003 Marseille 

and 
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National Office of Aerospace Studies and Research 
29, avenue de la Division Leclerc 
92320 Chatillon Sous Bagneux. 
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Figure 1. Scheme of the flow 
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